Material limitations of carbon-nanotube inertial balances: on the possibility of 
intrinsic yoctogram mass resolution at room temperature 
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We present a theoretical study of the intrinsic quality factor of the fundamental flexural vibration 
in a carbon nanotube and its dependence on temperature, radius, length and tension. In particular, 
we examine three- and four-phonon decays of the fundamental flexural mode within quantized elas- 
ticity theory. This analysis reveals design principles for the construction of ultrasensitive nanotube 
mass sensors: under tensions close to the elastic limit, intrinsic losses allow for single yoctogram 
mass resolution at room temperature, while cooling opens the possibility of sub-yoctogram mass 
resolution. 
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Recently, much experimental effort has gone into fab- 
ricating nano-electromechanical systems (NEMS) which 
employ carbon nanotubes as mechanical resonators [U-0] . 
The combination of small mass density and large mechan- 
ical stiffness makes carbon nanotube NEMS ideal candi- 
dates for ultrasensitive mass detectors which do not re- 
quire disruptive ionization of the analyzed molecule: re- 
cent experiments achieved a mass resolution of ~ 100 yg 
(1 yg=10~ 24 g)[H, However, for isotopic or chemi- 
cal identification of molecules adsorbed on the nanotube 
or real-time monitoring of chemical reactions yoctogram 
mass resolution is necessary. 

The mass resolution of nanotube NEMS depends 
strongly on the quality factor of the lowest flexural mode, 
which imposes a lower bound on the frequency differences 
which can be resolved. The quality factors found in most 
experiments have been surprisingly low, not exceeding 
Q fa 2000 |1 H0. Only very recently, Hiittel et al. 
measured Q > 10 5 , but at milikelvin temperatures. 

The intrinsic quality factor of a nanotube sets aside 
all extrinsic losses, such as defect or clamping losses, and 
sets an upper limit to the quality factor achievable in ex- 
periment. Both molecular dy namics simulations Q and 
analytical approaches 0, [l(| have been used to study 
the intrinsic quality factor theoretically, finding relatively 
low values similar to those in the experiments. However, 
simulations arc limited to extremely short tubes, and 
the analytical approaches to date have either employed 
simplified phonon spectra^ or focused solely on low- 
temperature four-phonon decays of the flexural mode [13] 
using an analysis appropriate to extremely long tubes. 

In this Letter, we present a theoretical analysis appro- 
priate to the length of tubes and operating temperatures 
commonly found in experiments. We discover that appli- 
cation of tension drastically increases the intrinsic qual- 
ity factor of nanotube oscillators and opens for the first 
time the theoretical possibility of single yoctogram mass 
resolution in such oscillators at room temperature, with 
further improvements possible with cooling. Indeed, the 
recent experiments of Wei et al. demonstrate a promising 
technique for controlling the tension in nanotubes (Tl|. 



For our analysis, we employ continuum elastic theory, 
which reliably describes long wavelength phonons in nan- 
otubes, to study the decay of the lowest flexural mode 
due to phonon-phonon interactions, the most important 
source of intrinsic losses in semiconducting, and possibly 
also metallic, nanotubes [lOj. Following the work of De 
Martino et al. [10] and Suzuura et al. [12j, we describe a 
nanotube as a rolled-up two dimensional elastic sheet and 
expand the free energy in powers of the strain tensor «y 
and the curvature tensor Sy . Curvature contributions to 
the free energy are generally small, but need to be in- 
cluded in the quadratic part to describe optical phonons 
properly. Using the isotropy of the elastic sheet leads to 



U 3 = Ci(Tru) 3 + C2Trw 2 Tra, 



U 2 = -(Tru) 2 + fiTru 2 



(Tr(S-So)Y 



(1) 
(2) 



where Ui and denote contributions to the energy den- 
sity including quadratic and cubic powers of and Sij . 
Also, A and fi denote the Lame moduli, k c the bending 
stiffness and Sq the curvature tensor of the undeformed 
tube, (i and (2 are anharmonic elastic constants. 

Next, using the standard expressions for iiij and Sij 
in cylindrical coordinates [10 . 1 1 21 ] , we express the poten- 
tial energy in terms of displacements u n with n = x,y, z 
(x denotes the tangential direction, y is along the tube 



axis and z is the radial direction) 



contains both lin- 



ear and nonlinear contributions in u n , which gives rise 
to two types of nonlinearities in the resulting Hamilto- 
nian: material nonlinearities, described by the explic- 
itly higher-order terms in Eq. ([5} , and geometric nonlin- 
earities, where the lower order terms, such as terms in 
Eq. ([1]), contribute at higher orders due to higher order 
terms in the expansion of in terms of the degrees of 
freedom u n . Material nonlinearities are often suppressed 
by a factor of kR (k being the phonon wave vector and 
R the tube radius) compared to geometric nonlinearities 
and can therefore be neglected for long wavelength 
phonons. However, when evaluating the three-phonon 
contribution to the quality factor, high energy phonons, 
whose wavelength is no longer comparable to the length 
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A [J/m 2 ] 


A* [J/m 2 ] 


[J/m 2 ] 


C2 [J/m 2 ] 


fc c [eV] 


59.57 


147.94 


-145.17 


-387.93 


1.46 



TABLE I: Lame moduli, anharmonic elastic constants of U3 
and bending stiffness for semiconducting carbon nanotubes 
from ab initio density-functional calculations. 
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FIG. 1: Energy of phonons versus kR for a R — 0.5 nm tube. 
The black dots denote phonon crossings, where three-phonon 
decays of the fundamental flexural mode are allowed. 



of the tube, play a crucial role. We, therefore, also con- 
sider in this work the contribution of material nonlinear- 
ities in the three-phonon decay channel. 

To determine accurate numerical values of the elastic 
constants in Eqs. ([1]) and we evaluate total energies 
of strained graphene sheets and various semiconducting 
nanotubes within density-functional theory in the local 
density approximation [13| and fitted the results to the 
continuum theory described above. Table U summarizes 
our results, which for the elastic constants of U2 are in 
good agreement with previous calculations [3, EH . 

We next impose canonical commutation relations to 
quantize the resulting elastic theory and diagonalize the 
quadratic part of the Hamiltonian by expressing the dis- 
placement operator in terms of normal modes. 

Figure [T] shows the resulting phonon dispersion curves 
for a tube with radius R = 0.5 nm, which are in good 
agreement with force constant models [l6[ and ab initio 
calculations (l5j |. In particular, we find as the lowest- 
frequency long wavelength modes two degenerate flexural 
modes with angular momentum component I = ±1 along 
the tube axis and a quadratic dispersion relation 



0{k 4 



n p M (\ + 2fj,) 



(3) 



where ujp denotes the angular frequency of the flexural 
phonon and pu = 7.68 x 10~ 7 kg/m 2 is the areal mass 
density of graphene. We also find a variety of low-lying 



optical phonons. Note that inclusion of curvature terms 
into Eq. ([I]) is crucial to obtain the correct optical phonon 
spectrum. In particular, the gap of the lowest optical 
branch vanishes if curvature terms are neglected. 

Next, we employ conservation laws to analyze the 
possible decay mechanisms of the fundamental flexural 
mode. In addition to energy conservation, translational 
invariance along the tube axis imposes conservation of k. 
Rotational invariance around the tube axis imposes con- 
servation of I. These conservation laws forbid the decay 
of the fundamental flexural mode in a three-phonon pro- 
cess which involves two other low- frequency modes Toj j . 
If, however, the fundamental flexural mode, which has 
a wavelength comparable to the length of the tube, in- 
teracts with a high-energy phonon of momentum k^ , en- 
ergy and momentum conservation enforce that the third 
phonon also has a high energy and momentum fc„ w k^. 
Because the flexural mode carries t = ±1, the angular 
momenta of the two high-energy phonons must differ by 
one. In sum, three-phonon decays of the fundamental 
flexural mode are only possible at crossings of two phonon 
bands whose angular momentum quantum numbers dif- 
fer by one. Inspection of Fig. Q] reveals that very few 
such crossings for energies comparable and smaller than 
ksT w 24 meV (at room temperature) exist. 

To compute the three-phonon contribution to the qual- 
ity factor, we evaluate the imaginary part of the Mat- 
subara Green function obtained from the lowest order 
bubble diagram, which contains a sum over intermediate 
high-energy phonon momenta. Because the energy un- 
certainty resulting from the short lifetimes of these 
modes is larger than the energy difference of phonons 
at neighboring allowed wave vectors of the finite-length 
tube, the sum can be converted into an integral. Then 
we use the energy-conserving <5-function to reduce the 
integral into a sum over allowed crossings. Our final ex- 
pression for the inverse quality factor resulting from the 
decay of a long wavelength flexural mode of wave vector 
k in a three-phonon process involving two high-energy 
modes (labeled p and v) at a phonon crossing at wave 
vector fc x and angular frequency u; x is 



^ = ^29F^(k,k x . 



/3h 2 n(u x )[n(uj x ) + 
RuF(k)uj^ \v£ - vl 



(4) 



where gF^v is a complicated coupling function that de- 
pends on the polarization vectors of all three phonons 
and gives a complete description of both material and 
geometric nonlinearities. Also, ^2 X denotes a sum over 
allowed crossings; v*, are the phonon group velocities 
and n{uS) — 1/ (exp(/3hui) — 1) is the Bose-Einstein factor 
with f3 = l/(fcsT) being the inverse thermal energy. 

Experiments with doubly clamped tubes inevitably in- 
volve some amount of strain. Here we consider the case 
of positive strain, corresponding to some amount of ten- 
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sion in the tube. To study the quality factor of such 
a strained tube, we expand the displacement around the 
new equilibrium value, taking into account the relaxation 
in the equilibrium radius. Evaluation of Eq. at the 
strained equilibrium configuration leads to two categories 
of additional contributions to the quadratic hamiltonian: 
one set of terms is already present in U2 and can be ab- 
sorbed into a redefinition of the linear elastic constants; 
the other set of terms gives the expected Hamiltonian for 
a string under tension, proportional to (d y u x ) 2 + (d y u z ) 2 
with a prefactor proportional to the tension. 

Computationally, having obtained the phonon disper- 
sions for a given radius, we numerically determine k x , 



and 



which are needed to evaluate gp^v and ul- 



timately the losses, for all relevant crossings. Our results 
indicate that geometric nonlinearities give the largest 
contribution to three-phonon losses, with material non- 
linearities contributing only about one percent. Also, we 
find that the primary effect from the application of ten- 
sion is to change the frequencies of the acoustic phonons, 
while the polarization vectors and the optical phonon fre- 
quencies change very little. 

Figure [2] shows our results for the three-phonon con- 
tribution to the inverse quality factors for tubes of typ- 
ical experimental radii as a function of inverse tempera- 
ture. At low temperatures [right side of Fig.[2Ja)], Q^ 1 
approaches zero exponentially, because the first Bosc- 
Einstein factor in Eq. (Q} rapidly diminishes the occu- 
pation of the high-energy modes at the crossings, which 
are the modes responsible for the scattering. At higher 
temperatures, T > 100 K, the modes associated with the 
relevant crossings are classically occupied and we find 
that Q^ 1 is proportional to temperature. 

Comparing the inverse quality factors of tubes of differ- 
ent radii in Fig.[2ja), we find that at low temperatures the 
tube with the smallest radius, R = 0.3 nm, exhibits the 
lowest dissipation, while at high temperatures its losses 
are largest [see Fig. [2fb)]. 

To understand this nontrivial radius dependence of Q3, 
we note that for T < 40 K the largest contribution to the 
losses comes from the crossing which is lowest in energy. 
The energy of this crossing of the flexural mode and the 
lowest optical mode is approximately equal to the gap of 
the optical mode given by E gap — y/k c / (5pAi)Qfi/R 2 , 
which depends sensitively on the radius of the tube. 
Therefore, at low temperatures, the occupation of the 
high-energy modes at this crossing is much more strongly 
suppressed for tubes with smaller radii (and therefore 
higher optical frequency), which leads to smaller losses 
according to Eq. (HJ). At higher temperatures, the con- 
tributions from other crossings become important (see 
Fig. [TJ. We find that those contributions depend quite 
sensitively on _R, which leads to the observed crossover 
behavior of the R = 0.3 nm tube. 

Figure OUa) shows the dependence of Q3 on the length 
L of the nanotube. In a tensionless tube (solid line), 
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FIG. 2: Log-log plots of the temperature dependence of Q3 
for tubes of different radii and length L = 500 nm. 
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FIG. 3: (a) Length dependence of Q3 for tubes with R = 
0.5 nm at T — 300 K. (b) Length dependence of mass resolu- 
tion 8m for tubes with R = 0.5 nm at T — 300 K. 



we find a remarkable cancellation between the length de- 
pendence of the coupling function, gF^v oc k 2 for small k, 
and the length dependence of wj? cx k 2 in Eq. (H}, result- 
ing in a quality factor which is insensitive to tube length 
beyond ~ 300 nm. Strained tubes do not exhibit this 
cancellation because ujp(k) is shifted by a constant pro- 
portional to the tension r if the tension is small (dashed 
line) (l8j . If r is large, ujp(k) oc \frk and Q3 becomes 
linear in L for long tubes (dotted line) |l8[ . 

For tubes shorter than 300 nm, Fig. |3][a) shows a sig- 
nificant length dependence of Q3 even for tensionless 
tubes. In particular, in tubes with lengths of only a few 
nanometers Q3 is reduced by more than an order of mag- 
nitude. This may be related to the small quality factors, 
Q « 1500, found by Jiang et al. @, who model the decay 
of the fundamental flexural mode in a 3 nm long singly 
clamped tube via molecular dynamics simulations. 

Next, we compare the magnitude of the computed 
intrinsic quality factor to experimental findings. At 
T = 300 K, we find Q 3 > 5 x 10 4 , which is at least 



4 



one order of magnitude larger than experimental results, 
Qexp < 2000 [J 01, suggesting that it is worthwhile to 
continue improving the control of losses in experiments. 
We find that the resulting (intrinsic) mass resolution, 
5m = 2M/Q, of a tensionless nanotube mass sensor de- 
pends sensitively on the tube length L with a minimum 
of ~ 5 yg for very short tubes [Fig. [3jb)]. On the other 
hand, we find (i) that application of tension can reduce 
the mass resolution 5m to a single yoctoqram if a tension 
close to the elastic limit, t c s» 100 nN [l9(, is applied and 
(ii) that 5m becomes independent of L for long tubes, 
simplifying the design and fabrication of actual devices. 

At lower temperatures, as discussed above, three- 
phonon processes are exponentially suppressed due to 
the energy gap of the optical modes, and the resulting 
dissipation becomes much smaller than the experimental 
findings in the milikclvin range [3|. We, therefore, now 
move on to consider losses from four-phonon decays. 

To estimate the role of four-phonon processes, we com- 
pute the leading order contribution, the fishbone dia- 
gram, due to a quartic coupling between four low-energy 
flexural modes. Following De Martino et al. and our 
findings for the three-phonon case, we only take into ac- 
count quartic geometric nonlinearities resulting from re- 
placing Uij in Eq. fl} by its nonlinear part. The resulting 
expression for the four-phonon contribution to the qual- 
ity factor contains a triple sum over intermediate phonon 
momenta and is given by (20j 



n 1 



E E 



919293 fl,£2,C3=± 

n(£i^i)n(£ 2 iJ2M£3^3) 
ra(£iwi + £ 2 w 2 + 6w 3 ) ' 

where qi denote intermediate momenta, Ui = ujp(qi) and 
V denotes the coupling function. Here, rj is the inverse 
lifetime associated with the dressed phonon lines repre- 
senting the actual phonons in the system with loss. 

In evaluating Eq. (|5|), proper account of the finite 
length of the tube, which leads to a finite spacing of 
the allowed wave vectors, is of crucial importance. For 
long wavelength flexural phonons in tubes of experi- 
mental lengths, the energy uncertainty A.E n associated 
with the observed lifetimes is actually much smaller than 
the energy difference Ai?Afc of phonons at neighboring 
wave vectors: Expressing 77 in terms of the quality fac- 
tor, 77 = cjf/(2ttQ), we find AE v /AE A k = 1/(4ttQ), 
which is much smaller than unity at low temperatures 
Thus, the sums over intermediate momenta cannot 
be converted into integrals for the tubes in the experi- 
ments. Converting the sums into integrals, as De Mar- 
tino et al.[I3] do, is appropriate for much longer tubes, 
but leads to an underestimate of Q4 due to inclusion of 
processes which are not present in the experiments. 

To describe four-phonon decays in experimentally rel- 
evant nanotubes, we carry out numerically the discrete 



triple sum over intermediate momenta in Eq. ([5]) using 
the discrete frequencies of a finite-length doubly clamped 
beam and experimentally observed inverse lifetimes rj. 
The result then gives the contribution to the observed 
losses from four-phonon processes. The resulting con- 
tributions to the quality factor are Q4 = 6.6 x 10 s at 
T = 1 K and Q 4 = 1.6 x 10 14 at T = 0.01 K for a tube 
of length L = 800 nm and R = 1.5 ran [3]. This indi- 
cates that, at low temperatures, four-phonon decays give 
only a small contribution to the observed losses, which 
are Q sa 10 4 at T = 1 K and Q sa 10 5 at T = 0.01 K 
Q. Although small, these losses are still much greater 
than the three-phonon contributions: Q3 = 3.1 x 10 27 at 
1 K and even greater at 0.01 K. Our analysis, contrary to 
the aforementioned theoretical studies, suggests that the 
losses observed by Hiittel et al.Q are mostly extrinsic. 
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